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Abstract In this paper, we investigate the downlink multiple-input-multiple- 
output (MIMO) broadcast channels in which a base transceiver station (BTS) 
broadcasts multiple data streams to K MIMO mobile stations (MSs) simul¬ 
taneously. In order to maximize the weighted sum-rate (WSR) of the system 
subject to the transmitted power constraint, the design problem is to find 
the pre-coding matrices at BTS and the decoding matrices at MSs. However, 
such a design problem is typically a nonlinear and nonconvex optimization 
and, thus, it is quite hard to obtain the analytical solutions. To tackle with 
the mathematical difficulties, we propose an efficient stochastic optimization 
algorithm to optimize the transceiver matrices. Specifically, we utilize the lin¬ 
ear minimum mean square error (MMSE) Wiener filters at MSs. Then, we 
introduce the constrained particle swarm optimization (PSO) algorithm to 
jointly optimize the precoding and decoding matrices. Numerical experiments 
are exhibited to validate the effectiveness of the proposed algorithm in terms 
of convergence, computational complexity and total WSR. 
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1 Introduction 

Multiple-input multiple-output (MIMO) transmission has been adopted for 
a various wireless communications standards such as the long-term evolution 
(LTE) and WiMax since it enables the enhancement of communication re¬ 
liability and channel capacity in wireless networks without the requirement 
of additional bandwidth or power mm- The single user MIMO communica¬ 
tion channels have extensively investigated in the literature, see, e.g. [1] and 
references therein. Recently, the multiuser MIMO channels have drawn much 
research attention for efficiently utilizing the spectrum and improving the sys¬ 
tem capacity [5]. In this paper, we study an important scenario in which a base 
transceiver station (BTS) equipped with multiple antennas sends different in¬ 
formation streams to the multiple MSs, each equipped with multiple anten¬ 
nas. Such a system model is also well-known as downlink multi-user MIMO 
broadcast channels (BCs), in which the major challenge is how to suppress 
co-channel interference at MSs. 

To mitigate interference and enhance the transmission data rate in down¬ 
link multiuser MIMO wireless networks, several linear precoding strategies 
using different algorithms have been studied EHS]. It was shown in [miT] that 
dirty paper coding (DPC) can achieve the capacity region of the multiuser 
MIMO BCs. Although DPC is the optimal transmission scheme for MIMO 
BCs, it is high complex for practical implementation. The authors in [ZlIH] 
have also proved that block diagonalization (BD) schemes can achieve the 
suboptimal transmission with lower complexity. Basically, BD is a version of 
the zero-forcing method in which the precoders are designed to null inter-user 
interference. Hence, BD provides both high per-MS rates and high sum rates 
by combining the benefits of spatial multiplexing and multiuser MIMO sys¬ 
tems. However, to efficiently suppress interference, the BD schemes require 
the number of transmit antennas at BTS to be larger than the sum of receive 
antennas of all MSs M- This makes BD inapplicable to the systems where 
there are a large number of MIMO MSs or the deployment of the large num¬ 
ber of transmit antennas at BTS is not applicable. Alternative approaches 
are based on deterministic optimization algorithms to maximize the system 
channel capacity mm- Due to the nonconvex nature of the design problem, 
the globally optimal solutions are generally not guaranteed. In addition, the 
deterministic iterative algorithms highly rely on the step size and/or gradient 
computation. As opposed to the previously existing methods, in this paper, 
we directly solve the weighted sum-rate (WSR) maximization problem by em¬ 
ploying the stochastic global optimization approach namely particle swarm 
optimization (PSO). 

The PSO algorithm is known as a class of swarm intelligence methods for 
solving global optimization problems. Introduced by Kennedy and Eberhart 
in 1995 [IHj, the key idea of PSO was inspired by social behaviour of bird 
flocks and fish schools. PSO seeks the optimal solution by transferring in¬ 
formation between particles and retains the global search strategy based on 
swarms. PSO can be efficiently implemented and take a reasonably short time 
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to converge to the optimal solution m- Hence, PSO has been exploited in 
many technical applications [IMS]- Reference m applied PSO to hnd opti¬ 
mal beamforming for multiuser MIMO systems. The algorithm in [18] aims to 
improve the total signal-to-interference-plus-noise (SINK) at each receiver and 
restricts single data steam to be transmitted to each MS. Motivated by the 
previous works, we develop the PSO algorithm to search the optimal solution 
for the nonconvex and intractable optimization problems of linear transceiver 
designs in downlink multiuser MIMO interference channels. In particular, we 
first formulate the design problems of finding the linear precoding and decod¬ 
ing matrices which maximize the total WSR of the systems subject to the total 
power constraints. From the conventional PSO algorithm which is normally 
applied to unconstrained optimization problem, we developed the constrained 
PSO algorithm which keeps the search space on the feasible regions by pro¬ 
jecting the particles into feasible sets. The constrained PSO algorithm is used 
to find the percoding matrices at BTS while the decoding matrices at MSs are 
the linear minimum mean square error (MMSE) Wiener filter solutions. As 
compared to the BD method which imposes the constraint on the number of 
antennas, the proposed method can be applied for more general scenarios of 
downlink MIMO systems. The computational complexity of the proposed PSO 
algorithm is shown to be in polynomial-time. The performance of the proposed 
algorithm shall be validated by numerical simulation results in compared with 
the BD method which is proven to achieve the suboptimal transmission rnisj. 

The rest of this paper is organized as follows. Section [5] introduces the 
downlink multiuser MIMO interference channels and presents the use of BD 
method for the design of the precoders and decoding matrices. In Section |31 
we present the PSO algorithm to find the optimal transceivers. The numer¬ 
ical results are provided in Section S) Finally, Section [S] presents concluding 
remarks. 

Notation: Matrices and vectors are respectively denoted by boldface cap¬ 
ital and lowercase letters. The transposition and conjugate transposition of 
complex matrix X are and X^, respectively. The i-th column of matrix 
X is denoted by X{:,i). I and 0 are respectively identity and zero matri¬ 
ces with the appropriate dimensions. tr(.), rank(.) and E(.) indicate the trace, 
rank and expectation operators, respectively. | |x| I 2 denotes the Euclidean norm 
while ||A’||f presents the Frobenius norm, x ^ CJV{x,Rx) mean that x is a 
complex Gaussian random vector with means x and covariance R^- 


2 System Model and BD Method 

In this section, we first present the model of downlink MIMO interference 
channels. Then, we briefly introduce the BD scheme to handle interference 
and maximize WSR. 
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Fig. 1 A system model of downlink MIMO interference channels. 


2.1 System model for downlink MIMO interference channels 


A downlink MIMO interference wireless system is considered as in Fig. ??, in 
which one BTS equipped with Nt antennas sends signals to K mobile stations 
(MSs), each equipped with the Nr^. antennas, for k = 1,..., AT. For the sake 
of presentation and without loss of generality, we assume that all MSs are 
equipped with the same number of antennas Nr^. = Nr for all k. Each MS 
receives dk = d data streams. Let Xk G be the signal vector transmitted 

from BTS to the fc-th MS. Before the signal is transmitted over Nt antennas, 
a linear precoding matrix Fk & is applied on Xk and then takes sum 

K 

for all MSs to form a signal vector s G given by s = ^ FkXk- Hence, 

k=l 


the transmitted power at BTS can be defined by 


K 

P = E{||sf} = ^ tr }Ff) . 


( 1 ) 


Without loss of generality, we assume that EjaifeX^} = Id and, hence, the 
constraint on the transmitted power from Eq. (HD can be rewritten as 


K 

P = ^ tr (FfcPf) < Pmax (2) 

fc=l 

where Pmax is the maximum allowable transmitted power at the BTS. 

In this paper, we consider the static flat-fading MIMO channels where the 
channel matrix from the BTS to the fc-th MS is denoted hy Hk G 
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We assume that all users are co-located and, hence, large scale fading can be 
excluded m- The channel entries are independent and identical distributed 
(i.i.d.) complex Gaussian variables and their magnitudes follow the Rayleigh 
distribution. The channel is assumed to be a block-fading which remains un¬ 
changed for a frame duration and varies independently for every frame. Since 
all signals are broadcasted from BTS, each MS receives not only its desired 
signal but also unintended signals from the other MSs. The received signal at 
the fc-th MS can be expressed as 


Uk =HkS 


K 

^ ^ H]^FiX£ Hk 

desired signal ^ ^ noise 

inter-user interference 


(3) 


where yk G C^’’^ ^ is the received signal at the fc-th MS and Uk ~ CN (O, a'j.lNr) 
is a complex Gaussian noise vector at the A:-th MS. To recover the desired sig¬ 
nal, the MSs apply the decoding matrix Wk & i^N^xd received signal 

yielding 


Xk = WkVk 


K 


= W^HkFkXk+ Y. Wj^HkFeXi+W^Uk. 
' ■■ ^ 1=1,e^k ' r ^ 


desired signal 


inter-user interference 


(4) 


Channel capacity at the fc-th MS can be calculated as [20] 


Uk = log2 \ld + W^HkFkF^Hj^WkR-,^ 


(5) 


where —W^ 


f K 

E HkFfF 

.i=\,t^k 


HuH < 

fe + Nr 


Wk is the correlation ma¬ 


trix of interference and noise in Eq.Q. The total WSR over K MSs can be 
expressed as 


K 

n = '^ujkUk ( 6 ) 

k^l 

where Wfc > 0 is the weight factor to present the different priority for the 
MSs. The problem of interest is to design the precoding matrices Fk and 
post-processing matrices Wk vci order to obtain the maximum WSR over K 
MSs. 


2.2 Block diagonalization for downlink multiuser MIMO interference channels 

An elRcient method to cancel the interference in MIMO downlink BCs is block 
diagonalization [3 El- This subsection introduces the block diagonalization 
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method for maximizing WSR which will be used as a benchmark for perfor¬ 
mance comparison. The key idea of the BD strategy is to find the precoding 
matrices and decoding matrices to force all inter-user interference to 
zero and maximize the total channel capacity. From Eq. o, the interference 
rejection condition must be expressed as 

HkFi=Q, V/c,andV£^fc. (7) 

Let us define the interference matrix as follows 

,Hl]^ (8) 

where “Hk G Then, Eq. ([7]) can be equivalently rewritten as 

UkFk = 0, Vfc. (9) 

It means that Fk has to lie in the null space oi'Hk- For the existence oi Fk 
satisfying Eq. we must have Nt — {K — l)Nr > d. The singular value 
decomposition (SVD) of 'Hk can be expressed as 

Uk = UkSkV^ ( 10 ) 

where Fk is the diagonal matrix with diagonal elements being singular values 
in decreasing order. Define Bk & as the matrix of the last d columns in 

matrix Vk, i-c., Bk lies on the null space of Bk- Thus, the precoder with Bk 
can make each block HkBk of the k-th MS free from inter-user interference. 
In order to separate the block channel HkBk into Nj. parallel sub-channels, 
an SVD is applied to HkBk as follows 

HkBk = UkAkV^ ( 11 ) 

where Ak is the diagonal matrix of decreasing singular values of HkBk- Define 
Dk £ as a matrix of the first d columns of Vk- Now, the precoder is 

defined as 

Fk=BkDkPl, (12) 

where Pk = diag[pk,i,Pk, 2 , ■ ■ ■ ,Pk,d] is the power allocation matrix. The de¬ 
coding matrix IV/c is chosen as the first d columns oiJJk- Eq. becomes 

Xk=W^HkBkDkPlxk+W^nk- (13) 

Then, channel capacity at the /c-th MS using BD can be expressed as 

P-kMD = y^log2(l H - ^^2 Pk,i)- (14) 

• T ^ h 

1—1 ^ 

The total WSR over K MSs can be expressed as 

K K d .2.. .N 

H ^ ^ log2(l H-^2^— Pk,i) 

k—1 k—1i—1 ^ 


( 15 ) 
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In order to maximize the WSR in Eq. (Ha, we need to find the optimal Pk by 
solving the following optimization problem 

K d f' 

max V] ^ w/c log 2 (1 + 2 ^ Pfc,i) (16a) 

K d 

-t. EE Pk,i ^ ^max- (16b) 

k—1i—1 

The Lagrangian of problem (fTBl) can be obtained as 

^Pk,i,>^) = - E E ^°S2(l + ^^^2 (17) 

Ai—1 "i—1 ^ 

K d 

+ a(EE Pk,i Pmax^ 
k—1i—1 

where A is the Lagrangian dual variables associated with the power constraint 
at BTS. Since problem (HID is obviously convex optimization, the optimal 
solutions of problem (IT51) must satisfy the set of Karush-Kuhn-Tucker (KKT) 


—^k — -ZZ — 


+ )ln2 

K d 


+ A = 0 


Y.Y.Pkd=Pma. 

A > 0. 


k—1i—1 


(18a) 

(18b) 

(18c) 


Then, the solution is given by 


Pk,i = max 0, 


Wfc 


A In 2 Ak^(i,i) 


(19) 


Finally, the water-filling (WF) algorithm in [3] is applied to find optimal A 
under the total power constraint at BTS. 


3 PSO for Downlink Multiuser MIMO Interference Channels 

The purpose of the current paper is to directly maximize the WSR of the 
networks in ([SD by searching the optimal precoding and decoding matrices. In 
order to reduce the search space and complexity, we adopt the MMSE Wiener 
filter at MSs since they take the effects of noise into account and are typically 
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exploited at the receivers because of their simplicity unj. Using Eq. the 
decoding matrix Wk can be found by m as 

Wk = ^ H^F^. ( 20 ) 

Then, we find the precoding matrices F^hy solving the optimization problem 
as follows 


K 


k^l 
K 

S.t. y^tr(.Ffc-Ff) < Prr 


K 


k=l 


(21a) 

(21b) 


The problem in Eq. m is the NP-hard nonconvex optimization which renders 
the mathematical challenges in finding the optimal solutions by the determin¬ 
istic optimization. This motivates us to develop the PSO algorithm to search 
for the global solutions. 

The standard PSO is originally known as a stochastic global optimization 
for unconstrained optimization problems. The PSO algorithm follows the sim¬ 
ulation of social particle behaviors in which particles are coordinated to share 
information to help each particle moving in search space towards the global 
optimal solution. The algorithm is initialized with a population of random 
particles within the feasible region, and then each particle uses experience of 
its own best location and the social best position to adjust its trajectory to 
search for a globally optimal solution in each evolutionary step Denote 

^ = {Fk}t as a set of precoding matrices. At the begin of the algorithm, 
the S particles are randomly initialized. The position and velocity of the i-th 
particle can be represented by Xi and Vi, respectively. The fitness of each 
partible can be calculated according to the cost function of the optimization 
problem. For each particle of the warm, the best previous moved position of 
the particle i is denoted as its best individual position Xpi,est,i- The best posi¬ 
tion of the entire swarm can be represented by Xgf,est ■ Based on the cognitive 
and social information, at the K-th iteration, the velocity and position of the 
particle i is adapted according to two following equations [ISlIig: 


V 


(k+1) 




(^) ( vM 


>'{k) \ 




y^gbest b 

^(k+ 1) _ _j_ ■^('^+1) 


(«)\ 


( 22 ) 

(23) 


where cq is called the inertial weight that controls the effect of the previous 
velocity of a particle on its current one. ci and C 2 are known as the cognitive 
and social parameters, respectively, which control the maximum step size. The 
random numbers ri^i and r 2 ,i are independently uniformly distributed in the 
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interval [0,1]. It is critical to choose these parameters since they affect the 
convergence characteristic of the PSO algorithm. An early convergence may 
occur with too small values of cq while a slow convergence may be result from 
Co being too large. In this paper, we select the values of cq = 0.7, ci = 1.494 
and C 2 = 1.494 according to experimental tests in [53] where these chosen 
parameters can commonly offer the good convergence characteristic. 

Although Eq. (1551) and (1531) perform a global search for unconstrained opti¬ 
mization problem, the design problem m is the constrained optimization one. 
Therefore, it is essential to develop a new constrained PSO algorithm. To deal 
with constraints, there are typically two approaches, namely the projection and 
penalty methods [53]. The basic idea of projection method is to project the 
solution of the associated unconstrained optimization problem into the feasi¬ 
ble region of the constraints while the penalty method focuses on transforming 
the constrained optimization problem into unconstrained optimization one by 
using the penalty function. Since choosing appropriate penalty parameters is 
difficult, this paper integrates the projection methods into the standard un¬ 
constrained PSO algorithm for finding the optimal precoders. We define the 
feasible region T of the optimization problem m as 

T=i^X = . (24) 

If particle Xi is not in F, then we project it back into the feasible region F. 
The projection of Xi into F is defined as 

X, = n{Xi) (25) 

which is mathematically expressed as 

K 

-ff({i^fc}f=i) = arg min (26) 

where Xi = {Fk}L ^ is the closed point in to Eq. (1531) is equivalently 
rewritten as 

K 

min '^tr{{Fk - Fk){Fk - Fk)^) (27a) 

K 

s.t. y^^trjFkFk) < Pmax- (27b) 

It is obvious that problem ([57]) is convex optimization. The Lagrangian ex¬ 
pression of problem (l27l) is defined as 

K 

C{{h}LiF) =^tr((Ffc -Ffe)(Ffe -Ffe)«) 

K 

-I- /i(^ tr(FfcFj, ) - Pmax) 
k=l 


(28) 
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where /r is the Lagrangian dual variables associated with the power constraint. 
The optimal solutions of problem (l?7ll must satisfy the set of Karush-Kuhn- 
Tucker (KKT) conditions 

(l + ^)Ffc-Ffc=0, A: = 1,2,...,X 

K 

tT{FkFf, ) = Pmax 

k^l 

p > 0 

which immediately implies that 

Fk= -Ffc, k = l,2,...,K. (30) 

Consequently, the constrained PSO using the projection method for finding 
the optimal transceivers can be described in Algorithm [T] where ATmax is the 
maximum number of iterations. 


(29a) 

(29b) 

(29c) 


Algorithm 1 : PSO for Downlink MIMO Broadcast Channels 
1: Declare variables X = {Fk}^_^. 

2: Inputs: , <Tfc, -Pmax,-/^max, S, cq, ci, C2, where Xbd 

is the solution of BD method. 

3: Set the iteration index k = 0. 

4: Generate a swarm of random solutions satisfying the power constraints {Xi}f_ 2 , where 
S is the swarm size. Set Xi = Xbd- 
5: Initialize for each particle Xpifest,i = Pi = 0- 

6: Calculate Wk from Eq. II20II for k = 1,2, ..K. 

7: Evaluate the cost function of each particle {TZ 
8: while k < A'max do 
9: Select the leader Xgi,est- 

10: For each particle, update the velocity and it position by Eqs. and if^ . 

11: If Xi is not into the feasible region, then use Eq. (I30II to project Xi into the feasible 

region: Xi ^ n{Xi). 

12: Calculate from Eq. II20II for A: = 1,2, ..K. 

13: Evaluate the cost function of each particle from Eq. Il21all . 

14: Update the Xpf)est,i for each particle. 

15: Update k, = k. 1 then move to the next "Do" task 

16: end while 


It is worth noting from Algorithm [T] that the particle with the highest WSR 
is chosen at each iteration and, thus, the objective is non-decreasing over 
iteration, i.e., 

nxits^)>nixltst)- ( 31 ) 

However, initial particle positions affect the convergence speed of the PSO al¬ 
gorithm. In this paper, we propose the use of the solution of the BD method as 
an initial particle while other initial particles are randomly chosen. It means 












Title Suppressed Due to Excessive Length 


11 


that our algorithm aims at searching the precoders which offers the higher 
sum-rate than those of BD methods since the precoders designed by the BD 
method are suboptimal solutions [7]. It should be emphasized that a key ad¬ 
vantage of the proposed PSO algorithm over the BS schemes is that the PSO 
algorthim does not impose any restriction of the number of users and antennas. 
The speed of convergence of Algorithm [T] shall be illustrated in the simulation 
results. The computational complexity of Algorithm [T] mainly relies on the 
computation of the initial particle from the BD scheme and the evaluation of 
the objective function for a swarm of S particles. For the BD method, the 
computational complexity which consists of the SVD computation and water¬ 
filling algorithm is 0{KN^ + K'^cP) j^. The major computation complexity 
of the PSO algorithm consists of computing matrix multiplications, matrix in¬ 
version, and determinant in Eqs. (IS|) and (Infill . Note that computational com¬ 
plexity of computing the inversion and determinant oi N x N matrix is 0{N^). 
Thus, evaluating the objective function approximately yields the complexity 
of 0{K‘^Nf). With S particles and the maximum number of iterations ATmax, 
the major computational complexity of Algorithm [T] is 0{SK^sy^K^N^). 


4 Simulation Results 

In this section, numerical simulation results are provided to evaluate the per¬ 
formance of our proposed PSO algorithm for downlink multiuser MIMO inter¬ 
ference channels. The performance of the proposed method is compared to the 
BD method |F| in terms of the total WSR. We denote the K user MIMO down¬ 
link channels with Nt transmit antennas at the BTS and W receive antennas 
at each MS by (W x Nr,d)^ |1S]. In simulations, noise variances are normal¬ 
ized = cr^ = 1. The Rayleigh fading channel coefficients are generated from 
the complex Gaussian distribution CA/"(0,1). Thus, we define SNR = All 
the numerical results are averaged over the 200 channel realizations. In the fol¬ 
lowing experiments, we select empirically the number of swarm size S = 700, 
hence, the local global optimum can be guaranteed and if S is chosen larger 
or S' —>■ oo, we can obtain the global solution. 

First, we investigate the convergence characteristic of the proposed PSO 
algorithm. We run the simulation for a (6 x 2,1)^ system for a random chan¬ 
nel realization. Weighed factors for 3 users are [wi, 072 ,^ 3 ] = [0.1, 0.2, 0.7]. The 
evolution of the objective function (weighted sum-rate) over iterations is illus¬ 
trated in Fig. As can been seen from Fig. the total WSR is significantly 
improved after several first iterations. The algorithm is converged to a fixed 
point in less than 150 iterations. 

Now, we compare the WSR of our proposed algorithm with that of the 
BD method. First, we consider a (6 x 2, {1,2})^ system where the condition 
Nt — {K — l)Nr > d is satisfied. Thus, the BD scheme can completely cancel 
inter-user interference. The WSR is plotted in Fig.|31 It has been revealed from 
Fig.[3]that the proposed algorithm offers a WSR performance improvement as 
compared to the BD method. The sum rate improvement is significant at low 
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Fig. 2 The convergence behavior of the proposed algorithm. 



Fig. 3 The average weighted sum-rate versus SNR for (6 X 2, {1, 2})® systems. 


SNR while the performance gap is smaller at high SNR. The reason is that at 
high SNR the system is interference-limited and the BD scheme approaches 
the optimal solution since it can cancel all interference. 

We also evaluate the WSR performance for different systems of (4 x 2, {1, 2})^. 
The average WSR is shown in Fig. U) It can be observed that the WSR per¬ 
formance of the BD method is not increased with an increasing SNR at the 
interference-limited region since the BD method cannot completely cancel all 
interference for these systems with Nt — {K — 1)W < d. In contrast, our 
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proposed algorithm directly maximizes the WSR and does not impose the re¬ 
striction on the number of antennas and MSs. It can been seen from Fig. 0] 
that our method consistently outperforms the BD method, especially for high 
SNR for both cases of c? = 1 and d = 2. 



Fig. 4 The average weighted sum-rate versus SNR for (4 X 2, {1, 2})® systems. 


5 Conclusion 

In this paper, we have studied the joint optimal transceiver designs for down¬ 
link multiuser MIMO interference channels. We have formulated the problem 
of finding the precoding and decoding matrices to achieve the maximization of 
the weighted sum rate subject to the total power constraints imposed on the 
BTS. The proposed approach is to employ an efficient stochastic PSO algo¬ 
rithm to obtain the optimal solutions. The numerical results indicate that the 
proposed PSO algorithm outperforms the BD method in terms of weighted 
sum rate. The weighted sum rate performance improvement is highly signif¬ 
icant when the number of users or receive antennas increases and the BD 
method cannot completely cancel interference. In addition, the numerical ex¬ 
periments have shown that the proposed methods can converge to a fixed point 
in less than a hundred and fifty iterations while the computational complexity 
per iteration is relatively low and affordable in practical applications. 
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